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Mapa míst konferencí ROBUST 


ONDŘEJ VENCÁLEK (CZ) a PAVEL STŘÍŽ (CZ) 


Abstrakt. Článek je úvodem do práce s mapami s dodatečnými popisky ve výpočetním 
prostředí R. Během naší práce používáme data z Eurostatu plus interní informace o kon- 
ferencích ROBUST. Vytvořená mapa je zobrazena v R i přes Rscript, vyexportována do 
png a pdf určeném pro Informační bulletin na tisk, navíc zobrazena ve webovém prohlížeči 
offline i online za pomocí knihovny plotly a jejich prostředí Chart Studio určeném pro práci 
v týmech a zaměřeném na zveřejňování grafů. 


Klíčová slova. programovací jazyk R, Eurostat, sf, dplyr, ggplot2, ggrepel, plotly, Chart 
Studio, teltk. 


A MAP OF ROBUST CONFERENCE VENUES 


Abstract. The article is a brief introduction to a map creation in R programming lan- 
guage with additional information. We primarily use data from the Eurostat servers and 
internal information about the ROBUST conference series. A preview of the map is shown 
in R and Rscript environments, exported to the png image and the pdf file for printing 
purposes of the Information Bulletin as well as previewed offline and online using the plotly 
library and their Chart Studio designed for collaborative work and chart publishing. 


Keywords. R programming language, Eurostat, sf, dplyr, ggplot2, ggrepel, plotly, Chart Stu- 
dio, teltk. 


1. Úvod 


Konferenční řada ROBUST má v roce 2020 20+20. výročí, a proto jsme se roz- 
hodli vytvořit tzv. „robust(n)í mapu“, tj. mapu Česko-Slovenska se zvýrazněním 
míst. V neformálním duchu jako jsou tyto konference, možná jako dárek s cílem 
si zkusit pár fint v R. 


2. Data z Eurostatu 


Z webové stránky https://ec.europa.eu/eurostat/web/gisco/geodata/re- 
ference-data/administrative-units-statistical-units/nuts stáhneme si 


soubor zip v řádku NUTS 2016, 14/03/2019, 1:1 Million, sloupec SHP. NUTS ne- 
boli Nomenklatura územních statistických jednotek jsou územní celky vytvořené 
pro statistické účely Eurostatu. NUTS 0 odpovídá státu, NUTS 1 území (ČR), 
NUTS 2 regionům, NUTS 3 krajům, LAU 1 (dříve NUTS 4) okresům a LAU 2 
(dříve NUTS 5) obcím. 
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Užijeme adresář, kde budeme mít zdrojový kód v R. Konkrétně: 


$ wget https://ec.europa.eu/eurostat/cache/GISCO/distribution/N 
> v2/nuts/download/ref-nuts-2016-01m.shp.zip 


Ze staženého souboru potřebujeme extrahovat jen jediný. Na vykreslení mapy 
nám stačí. 


$ unzip ref-nuts-2016-01m.shp.zip NUTS RG 01M 2016 4326.shp.zip 


RG značí regiony (angl. regions; vedle toho najdeme BN/boundaries pro hra- 
nice a LB/labels pro popisky), 01M je vybrané měřítko, 2016 je rok regulace sta- 
tistické jednotky (aktualizace 14.3.2019) a 4326 značí EPSG:4326, také známé 
jako geografická projekce WGS84 v desetinách stupňů. LEVL odpovídá NUTS. 

Ostatní soubory jsou v dalších geografických projekcích (EPSG:3035 pro Ev- 
ropskou projekci LAEA a EPSG:3857 pro WGS84 / Pseudo-Mercator), menších 
územních celcích a rozlišené dle regionů, hranic a popisků. 

Tento soubor kompletně rozbalíme: 


$ unzip NUTS RG 01M 2016 4326.shp.zip 


Získáme tak hlavní soubor vektorových dat (shp), indexový soubor (shx) a da- 
tabázový soubor (dbf). Tyto soubory jsou povinné, vedle toho získáme dva do- 
plňkové soubory: zaznačení kódové stránky pro dbf (cpg) a přesné informace 
o projekci (prj). 


3. Softwarová příprava 


Pokud nemáme, nainstalujeme si R a potřebné knihovny, v našem případě mini- 
málně libssl-dev, bude jej potřebovat plotly u interaktivní online verze: 


$ sudo apt install r-base-core 
$ sudo apt install libssl-dev 


Ve výpočetním prostředí R si doinstalujeme knihovny: 


$R 
> install.packages(c("sf","dplyr","ggplot2","ggrepel","plotly","tcltk")) 
> ga 


Ideální je instalovat knihovny jednu za druhou, skrz možné chybové zprávy. 
Dále si vytvoříme soubor robust-celek.R, který si postupně okomentujeme. 
Může se vám hodit TgXLive, ImageMagick nebo GraphicsMagick, prohlížeč 

rastrových obrázků (např. ristretto či GIMP) a prohlížeč pdf souborů (např. 

Okular či XpdfReader). Necháváme na laskavém čtenáři, jak hluboce se chce do 

článku ponořit. 
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4. Stati(sti)cká mapa 


Načteme si potřebné knihovny. 


1 library(sf); library(dplyr); library (ggplot2) 
2 library(ggrepel); library(tcltk); library(plotly) 


Připravíme si data o konferenční řadě v proměnné rob(ust). 


3 ft Data o ROBUSTech 4 

4 rob = data.frame( 

5 location = c("Načetín","Kost","Slavonice","Adršpach","Plasy", 

6 "Liblice" ,"Herbertov","Malenovice","Lednice","Radešín", 

7 "Nečtiny" ,"Hejnice","Třešť","Lhota nad Rohanovem","Pribylina", 

8 "Králíky" ,"Němčičky","Jetřichovice","Kurzovní","Rybník","Bardějov"), 
9 longitude = c(13.2634230,15.1350694,15.3195072,16.1025647,13.3838904, 
10 14.5853101,14.3358220,18.4153645,16.8023564,16.0841326, 

11 13.1450920,15.2190330,15.4885601,13.6832665,19.7801583, 

12 16.7519537,16.8280667 ,14.3837522,17.2096600,12.6724633,21.2713128), 
13 latitude = c(50.5489359,50.4902238,49.0149020,50.6203202,49.9340679, 
14 50.3237900,48.6200863,49.5678849,48.8082447 ,49.466635bb, 

15 49.9541085,50.8502057,49.3102300,49.1399375,49.1334503, 

16 50.1007519,48.9331628,50.8655207,50.0785817,49.5053812,49.2926572), 
17 Termín = c("letní","zimní","letní","zimní","letní", 

18 "zimní","letní","zimní","letní","zimní", 

19 "letní","zimní","letní","zimní","letní", 

20 "zimní","letní","zimní","letní","zimní","2020"), 

21 | roky = seg(1980,2020,by=2) 

22) 


Načteme si mapu a vytvoříme si základ v proměnné p(lotly). 


23 4 Základ mapy ff 

24 map = st read("NUTS RG 01M 2016 4326.shp", stringsAsFactors=FALSE) 
25 CZO = map %>% filter(NUTS ID %in/ c("CZ","SK")) 

26 p = ggplot(rob) + 

27. geom sf(data = CZO, show.legend = FALSE, color="black") + 
28. xlab("") + ylab("") + 

29 | geom point(data=rob, 

30 aes(x=longitude, y=latitude, fill=Termín), 

31 pch=21, size=4, alpha=I(0.7)) + 

32 scale fill manual(values=c("green","red","blue")) + 

33. theme void() + 


34 | theme( 

95 legend.position="bottom", 

36 legend.title = element blank(), 

37 panel.grid.major = element line(linetype="blank"), 
38 panel.grid.minor = element line(linetype="blank") 


39.) 
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Náš poslední krok je přidat popisky. 


Hit Přidání popisku, statická verze Hfft 
pfinal = p + geom label repel(data=rob, aes(x=longitude, y=latitude, 
label=location), cex=3) 


Nyní se můžeme podívat na výsledek: 


Hit Zobrazení v prostředí R fff 
pfinal 


Pokud bychom rádi mapu viděli přes Rscript, můžeme užít X11 () pro operační 
systém Linux, resp. windows() či guartz() pro Microsoft Windows a Mac OS 
X. Pozastavení zrealizujeme přes knihovnu tcltk. 


HH Zobrazení přes Rscript ff 

X11() £ windows() nebo guartz() 

pfinal 

prompt <- "Mezerníkem se zavře graf" 
capture <- tk messageBox(message = prompt) 


dev.off() 
Dále si ukážeme uložení do pdf. 
HX Uložení do PDF, verze tisková K 


cairo pdf ("Robust.pdf", width=10, height=7) 
pfinal 
dev.off() 


Pokud bychom si přáli ořezat ochrannou bílou zónu, lze to udělat přes příkaz 
$ pdfcrop --hires --margins 1 Robust.pdf 


Nástroj pdfcrop je součástí TEXové distribuce TpXLive. Nula se zpravidla 
nepoužívá kvůli přesahům Bézierových křivek v písmech. 

Jiná možnost je uložit si mapu do rastrového obrázku a ten si zobrazit přes 
prohlížeč obrázků. 


HH Uložení do PNG, verze na web fff 
png ("Robust.png") 

pfinal 

dev.off() 

browseURL ("Robust ..png") 


Pokud chceme ořezat ochrannou bílou zónu i zde, lze i to: 


$ gm convert Robust.png -flatten -fuzz 1/ -trim +repage Robust-crop.png 
$ £ convert Robust.png -flatten -fuzz 1/ -trim +repage Robust-crop.png 
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5. Interaktivní mapa offline 


Náš další krok je připravit si verzi pro internet. Inspiraci lze hledat na serveru 


https://plot.ly/r/ 


Prvně zasáhneme do vzhledu mapy. 


60 4+ Přidání popisků, interaktivní verze ff 

61 p = ggplot(rob) + 

62 | geom sf(data = CZO, show.legend = FALSE) + 

63. xlab("") + ylab("") + 

64 | geom point(data=rob, 

65 aes(x=longitude, y=latitude, fill=Termín, text=paste(roky,'<br>',location)), 
66 pch=21, size=4, alpha=I(0.7)) + 

67. | scale fill manual(values=c ("green","red","blue")) + 

68 | theme void() + 


69. theme( 

70 legend.position="bottom", 

71 panel.grid.major = element line(linetype="blank"), 
72 panel.grid.minor = element line(linetype="blank") 
73.) 


74 p <- ggplotly(p, tooltip= c("text")) £ c("location")) 
75 p <- p 44 

76. layout(legend = list( 

TT orientation = "h", x = 0.4, y =0 

78.) 


Zobrazit si mapu v R v offline režimu je jednoduché. 


79 4 Zobrazení mapy bez nutnosti účtu na plot.ly Hf 
80 p 


Na první dobrou se nám nepodařilo mapu zobrazit přes Rscript. Jeden ze 
zajímavých tipů na internetových fórech] bylo místo p použít jednu z možností: 
p 4>4 print 
print (p) 
show (p) 

To sice začalo fungovat, pokud jsme blok kódu volali v R přes příkaz source 
a v RStudiu, ale nikoliv přes Rscript. 

Rešení jsme objevili zde? Nemáme-li, doinstalujeme si webové prohlížeče: 
$ sudo apt install chromium-browser 
$ sudo apt install firefox 


Návod pro chromium-browser funguje bezvadně. 


INamátkou zde: https://stackoverflow.com/guestions/26643852/ggplot-plots-in- 
scripts-do-not-display-in-rstudio 

2https://stackoverflow.com/guestions/44048347/r-open-plotly-in- 
standalone-window 
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print app <- function(widget) 1 


temp <- paste(tempfile('plotly'), 'html', sep = '.') 
htmlwidgets::saveWidget (widget, temp, selfcontained = FALSE) 
system(sprintf ("chromium-browser -app=file:///s", temp)) 
* system(sprintf ("firefox file://4s", temp)) 
* system(sprintf ("sleep 5")) 
temp 

] 

print app(p) 


Ovšem u prohlížeče firefox se dočasný soubor smaže dřív, než jej prohlížeč 
dokáže otevřít. Rešení je pozastavit funkci print app a jsou to zakomentované 


řádky |85|a 


6. Interaktivní mapa online 


Náš poslední robustní úkol je mapa v online režimu. 


e Otevřeme sihttps://chart-studio.plot.ly/ 


Naklikneme Sign Up. 

Vyplníme potřebné registrační údaje. Dáme SIGN UP. 

Potvrdíme údaje přes odkaz v došlém emailu. 

Po najetí na vytvořené uživatelské jméno (vpravo nahoře) dáme Settings. 
Vlevo v menu zvolíme API Keys. 

Přes Generate Key si necháme vytvořit API klíč. 

Někam si uživatele a vytvořený API klíč uložme. Hned je použijeme. Heslo 
je dobré si též uložit, ale v R jej potřebovat nebudeme. 


Vrátíme se zpět do zdrojového kódu a zapišme: 


Hit Zobrazení mapy na serveru plot.ly ff 

Sys.setenv("plotly username"="VÁŠ UŽIVATEL") %£ <-- zasáhněte 
Sys.setenv("plotly api key" ="VÁŠ API KLÍČ") £ <-- zasáhněte 
chart link = api create(p, filename="Robust-2020") 

chart link 


Proběhlo-li vše v pořádku, otevře se webový prohlížeč a vznikly (nebo se ak- 
tualizovaly) dvě záložky: Robust-2020 Grid a Robust-2020. 

Pracovní verze mapy je dostupná přes stránku 
ly/create/?fid=vencalek:5%/, zkrácený odkaz je 

Nezávislý test je uložen na 
resp. 


Gratulujeme! Jsme hotoví. 
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7. Mapa na produkčním serveru 


Webové stránky ČStS — — jsou přepracované. Pokud z menu vy- 
bereme KONFERENCE a dále ROBUST, mapa se otevře v pravé horní části. 


Přímý odkaz na výsledek našich snah je http://www.statspol.cz/konference/ 
robust/. 


INFORMAČNÍ BULLETIN KONFERENCE KNIHOVNIČKA 


dt ROBUST 


Konference ROBUST je nejvýznamnějším 
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